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Abstract. Stack filters are a special case of non-linear filters. They have 
a good performance for filtering images with different types of noise while 
preserving edges and details. A stack filter decomposes an input image 
into several binary images according to a set of thresholds. Each binary 
image is then filtered by a Boolean function, which characterizes the 
filter. Adaptive stack filters can be designed to be optimal; they are 
computed from a pair of images consisting of an ideal noiseless image 
and its noisy version. In this work we study the performance of adaptive 
stack filters when they are applied to Synthetic Aperture Radar (SAR) 
images. This is done by evaluating the quality of the filtered images 
through the use of suitable image quality indexes and by measuring the 
classification accuracy of the resulting images. 
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1 Introduction 

SAR images are generated by a coherent illumination system and are affected 
by the coherent interference of the signal from the terrain [1 . This interference 
causes fluctuations of the detected intensity which varies from pixel to pixel, an 
effect called speckle noise, that also appears in ultrasound-B, laser and sonar 
imagery. Speckle noise, unlike noise in optical images, is neither Gaussian nor 
additive; it follows other distributions and is multiplicative. Classical techniques, 
therefore, lead to suboptimal results when applied to this kind of imagery. The 
physics of image formation leads to the following model: the observed data can 
be described by the random field Z ^ defined as the product of two independent 
random fields: X, the backscatter, and 1", the speckle noise. The backscatter is 
a physical magnitude that depends on the geometry and water content of the 
surface being imaged, as well as on the angle of incidence, frequency and po- 
larization of the electromagnetic radiation emitted by the radar. It is the main 
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source of information sought in SAR data. Different statistical distributions have 
been proposed in the literature for describing speckled data. In this work, since 
we are dealing with intensity format, we use the Gamma distribution, denoted 
by r, for the speckle, and the reciprocal of Gamma distribution, denoted by 
r^^, for the backscatter. These assumptions, and the independence between the 
fields, result in the intensity law for the return [2]. This family of distributions 
is indexed by three parameters: roughness a, scale 7, and the number of looks 
n, and it has been validated as an universal model for several types of targets. 
Speckle has a major impact on the accuracy of classification procedures, since 
it introduces a low signal-to-noise ratio. The effectiveness of techniques for com- 
bating speckle can be measured, among other quantities, through the accuracy 
of simple classification methods. The most widespread statistical classification 
technique is the Gaussian maximum likelihood classifator. Stack filters are a spe- 
cial case of non-linear filters. They have a good performance for filtering images 
with different types of noise while preserving edges and details. Some authors 
have studied these filters, and many methods have been developed for their con- 
struction and applicaton as in [5]. These filters decompose the input image, by 
thresholds, in binary slices. Each binary image is then filtered using a Boolean 
function evaluated on a sliding window. The resulting image is obtained sum- 
ming up all the filtered binary images. The main drawback in using stack filters 
is the need to compute optimal Boolean functions. Direct computation on the 
set of all Boolean functions is unfeasible, so most techniques rely on the use of a 
pair of images: the ideal and corrupted one. The functions are sought to provide 
the best estimator of the former using the latter as input. The stack filter design 
method used in this work is based on an algorithm proposed by Yoo et al. [4]. 

We study the application of this type of filter to SAR images, assessing its 
performance by evaluating the quality of the filtered images through the use of 
image quality indexes like the universal image quality index and the correlation 
measure index and by measuring the classification accuracy of the resulting 
images using maximum likelihood Gaussian classification. 

The structure of this paper is as follows: In Section [5] we summarise the Q° 
model for speckled data. Section |2] gives an introduction to stack filters, and 
describes the filter design method used in this work. In Section 2] we discuss 
the results of filtering through image quality assesment and classification perfor- 
mance. Finally, in Section [5] we present the conclusions. 

2 The Multiplicative Model 

Following [51, we will only present the univariate intensity case. Other formats 
(amplitude and complex) are treated in detail in [5]. 

The intensity Q'^ distribution that describes speckled return is characterized 
by the following density: 



L^r{L~a) z-^-i 
~ -f'^r{L)r{-a) (7-KLz)'f'-"' 
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where — a,7, z > 0, L > 1, denoted Q^{a,j,L). 

The a parameter corresponds to image roughness (or heterogenity). It adopts 
negative values, varying from — oo to 0. If a is near 0, then the image data are 
extremely heterogeneous (for example: urban areas), and if a is far from the 
origin then the data correspond to a homogeneous region (for example: pasture 
areas). The values for forests lay in-between. 

Many filters have been proposed in the literature for combating speckle noise, 
among them the ones by Lee and by Frost. These filters will be applied to speck- 
led data, along with the filter proposed in this work. For quality performance 
the comparision will be done between the stack filter and the Lee filter. Classi- 
fication performance will be assessed by classifying data filtered with the Lee, 
Frost and stack filters using a Gaussian maximum likelihood approach. 

3 Stack Filters 

This section is dedicated to a brief synthesis of stack filter definitions and design. 
For more details on this subject, see [41617) . 

Consider images of the form X : S {0, . . . ,M}, with S the support and 
{0, . . . , M} the set of admissible values. The threshold is the set of operators 
T™ : {0, . . . , M} -> {0, 1} given by 

^ ' \0 if X < m. 

We will use the notation X™ — T™ (x) . According to this definition, the value of 
a non-negative integer number x G {0, . . . , M} can be reconstructed making the 
summation of its thresholded values between and M . Let X = {xq, . . . , Xn-i) 
and Y = (j/o, • ■ • ,2/n-i) be binary vectors of length n, define an order relation 
given by X < y if and only if holds that Xi < yi for every i. This relation is 
refiexive, anti-symmetric and transitive, generating therefore a partial ordering 
on the set of binary vectors of fixed length. A boolean function /: {0,1}" — > 
{0, 1}, where n is the length of the input vectors, has the stacking property if 
and only if 

VX, Y e {0, 1}", X<Y^ f{X) < f{Y). 

We say that / is a positive boolean function if and only if it can be written 
by means of an expression that contains only non-complemented input variables. 
That is, f{xi,X2, ■ ■ ■ , Xn) — ViLi /\jePi -^i' where n is the number of arguments 
of the function, K is the number of terms of the expression and Pi is a subset of 
the interval {1, . . . , N}. 'Y' and '/\' are the AND and OR Boolean operators. It 
is possible to proof that this type of functions has the stacking property. 

A stack filter is defined by the function Sf. {0, . . . , M}" {0,...,M}, 
corresponding to the Positive Boolean function f{xi,X2,---,Xn) expressed in 
the given form by ([3]). The function Sf can be expressed by means of Sf{X) = 

E™=i/(r™(x)). 
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In this work we applied the stack filter generated with the fast algorithm 
described in [i]. 

Stack filters are built by a training process that generates a positive boolean 
function that preserves the stacking property. Originally, this training is per- 
formed providing two complete images on S, one degraded and one noiseless. 
The algorithm seeks the operator that best estimates the later using the former 
as input, and as a means of measuring error. 

The implementation developed for this work supports the application of the 
stack filter many times. Our approach consists of using a set of regions of interest, 
much smaller than the whole data set, and relying on the analysis the user makes 
of these information. Graphical and quantitative analyses are presented. The user 
is prompted with the mean value of each region as the default desired value, but 
he/she can choose other from a menu (including the median, the lower and 
upper quartiles and a free specification) . This freedom of choice is particularly 
useful when dealing with non-Gaussian degradation as is the case of, for instance, 
impulsive noise. 



4 Results 

In this section, we present the results of building stack filters by training. These 
filters are applied to both simulated and real data. The stack filters obtained are 
compared to SAR image filters. This comparision is done by assessing smooth- 
ing and edge preservation through image quality indexes and by evaluating the 
influence of filtering on classification performance. 



4.1 Image quality assesment 

The indexes used to evaluate the quality of the filtered images are the universal 
image quality index [S] and the correlation measure f3. The universal image 
quality index Q is given byequation ([1]) 



axY 2XY 2ax(JY 



(1) 



where - (N -l)-^Sl,{X, - X)^, ^ {N - l)-^St^{Y, - Y)\ X = 
N-^S^^^Xi and Y = N-^S^^^Y,. The dynamic range of index Q is [-1,1], 
being 1 the best value. To evaluate the index of the whole image, local indexes 
Qi are calculated for each pixel using a suitable square window, and then these 
results are averaged to yield the total image quality Q. The correlation measure 
is given by 

(2) 

where V^X and V^F are the Laplacians of images X and Y , respectively. 

In Table [T] the correlation measure and the quality index Q are shown. The 
comparison is made between Lee filtered and stack filtered SAR images. To this 
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Table 1. Statistics from image quality indexes 



/3 index Q index 





Stack filter 


Lee filter 


Stack filter 


Lee filter 


contrast 




^ sp 


Q SQ 


Q SQ 


10:1 
10:2 
10:4 
10:8 


0.1245 0.0156 
0.0964 0.0151 
0.0267 0.0119 
-0.0008 0.0099 


0.0833 0.0086 
0.0663 0.0079 
0.0421 0.0064 
0.0124 0.0064 


0.0159 0.0005 
0.0154 0.0005 
0.0124 0.0008 
0.0041 0.0013 


0.0156 0.0004 
0.0148 0.0004 
0.0120 0.0006 
0.0021 0.0006 



end, a Monte Carlo experiment was performed, generating 1000 independent 
replications of synthetic 1-look SAR images for each of four contrast ratios. The 
generated images consist of two regions separated by a vertical straight border. 
Each sample corresponds to a different contrast ratio, wich ranges from 10:1 
to 10:8. This was done in order to study the effect of the contrast ratio in the 
quality indexes considered. 

It can be seen that, according to the results obtained for the /? index, the 
stack filter exhibits a better performance at high contrast ratios, namely 10:1 
and 10:2, while the Lee filter shows the opposite behavior. The results for the Q 
index show slightly better results for the stack filter all over the range of contrast 
ratios. It is remarkable the small variance of these estimations, compared to the 
mean values obtained. 

Fig. [T] shows the boxplots of the observations summarized in Table [TJ From 
the plots for the /3 index, it can be seen that, the Lee filter has a lower degree of 
variability with contrast and that both are almost symmetric. The plots of the 
Q index show a better performance for the stack filter for all the contrast ratios 
considered. 



4.2 Classification performance 



The equality of the classification results are obtained by calculating the confusion 
matrix, after Gaussian Maximum Likelihood Classification (GMLC). 



Fig. 2(a) left, presents an image 128 x 128 pixels, simulated with two regions: 
samples from the 1.5, 7* ^ 5 1) and from the 10, 7* 1) laws form 
the left and right halves, respectively, where 7* „ denotes the scale parameter 
that, for a given roghness a and number of looks n yields an unitary mean 



law. In this manner. Fig. 2(a) presents data that are hard to classify: extremely 
heterogeneous and homogeneous areas with the same mean, with the lowest 
possible signal-to- noise ratio (n = 1). The mean value of the dashed area was 



used as the "ideal" image. Fig. 2(b) and 2(c) left, show the result of applying 
the resulting filter once and 95 times, respectively. The right side of Fig. 2(a) 



Fig. 2(b) and Fig. 2(c) present the GMLC of each image. Not only the pointwise 



improvement is notorious, but the edge presevation is also noteworthy, specially 



in Fig. 2(c) right, where the straight border has been completely retrieved. 
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(a) Values of /3, Lee filter 




(b) Values of /3, Stack filter 




(c) Values of Q, Lee filter (d) Values of Q, Stack filter 

Fig. 1. Boxplots of the quality indexes 



Fig. 13] compares the performance of the proposed stack filter with respect to 



two widely used SAR filters: Lee and Kuan. Fig. 3(a) presents the original data, 
and the regions of interest used for estimating the Boolean function. In this case, 
again, the mean on each region was used as the 'ideal' image. Fig. |3(b)| Fig 3(c) 



and Fig. [3(d)] present the result of applying the Frost, Lee and Stack filters (one 
and 22 iterations) to the original SAR data. The right side of previous figures 
present the corresponding GMLC. The stack filter produces better results than 
classical despeckling techniques. 

Table [5] presents the main results from the confusion matrices of all the 
GMLC, including the results presented in [9] which used the classical stack filter 
estimation with whole images. It shows the percentage of pixels that was labeled 
by the user as from region i?,; that was correctly classified as belonging to region 

for 1 < I < 3. "None" denotes the results on the original, unflltered, data, 
"Sample Stack /c" denotes our proposal of building stack filters with samples, 
applied k times, "Stack fc" the classical construction applied k times, and "Frost" 
and "Lee" the classical speckle reduction filters. 
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(a) Simulated image and (b) One iteration and (c) 95 iterations and 
GMLC GMLC GMLC 



Fig. 2. Training by region of interest: simulated data 




(a) Image, samples and (b) Frost and GMLC (c) Lee and GMLC 

GMLC 




(d) Stack Filter 20 and 
GMLC 



Fig. 3. Training by region of interest: real image 

It is clear the superior performance of stack filters (both classical and by 
training) over speckle filters, though the stack filter by training requires more 
than a single iteration to outperform the last ones. 

Stack filters by training require about two orders of time less than classical 
stack filters to be built, and they produce comparable results. Using regions of 
interest is, therefore, a competitive approach. 

5 Conclusions 

In this work, the effect of adaptive stack filtering on SAR images was assessed. 
Two viewpoints were considered: a classification performance viewpoint and a 
quality perception viewpoint. For the first approach, the Frost and Lee filters 
were compared with the iterated stack filter using a metric extracted from the 
confusion matrix. A real SAR image was used in this case. For the second ap- 
proach, a Monte Carlo experience was carried out in which 1-look synthetic SAR, 
i.e., the noisiest images. w(;re generatcxl. In this case, the Lee filter and a one 
pass stack filter were compared for various degrees of contrast. The (3 and the Q 
indexes were used as measures of perceptual quality. The results of the /3 index 
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Table 2. Statistics from the confusion matrices 



Filter 




R2/R2 


Ri/Rz 


None 


13.40 


48.16 


88.90 


Sample Stack 1 
Sample Stack 22 


9.38 
63.52 


65.00 
74.87 


93.19 
96.5 


Stack 1 
Stack 40 
Stack 95 


14.35 
62.81 
63.01 


64.65 
89.09 
93.20 


90.86 
94.11 
94.04 


Frost 
Lee 


16.55 
16.38 


55.54 
52.72 


90.17 
89.21 



shows that the stack filter performs better in cases of high contrast. The results 
of the Q index show slightly better performance of the stack filter over the Lee 
filter. This quality assessment is not conclusive but indicates the potential of 
stack filters in SAR image processing for visual analysis. The classification re- 
sults and the quality perception results suggest that stack filters are promising 
tools in SAR image processing and analysis. 
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